A solar cycle lost in 1793-1800: 
Early sunspot observations resolve the old mystery 
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ABSTRACT 

Because of the lack of reliable sunspot observation, the quality of sunspot 
number series is poor in the late 18th century, leading to the abnormally long 
solar cycle (1784-1799) before the Dalton minimum. Using the newly recovered 
solar drawings by the 18-19th century observers Staudacher and Hamilton, we 
construct the solar butterfly diagram, i.e. the latitudinal distribution of sunspots 
in the 1790's. The sudden, systematic occurrence of sunspots at high solar lati- 
tudes in 1793-1796 unambiguously shows that a new cycle started in 1793, which 
was lost in traditional Wolf's sunspot series. This finally confirms the existence 
of the lost cycle that has been proposed earlier, thus resolving an old mystery. 
This letter brings the attention of the scientific community to the need of revising 
the sunspot series in the 18th century. The presence of a new short, asymmetric 
cycle implies changes and constraints to sunspot cycle statistics, solar activity 
predictions, solar dynamo theories as well as for solar-terrestrial relations. 

Subject headings: Sun: activity — sunspots 
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Introduction 



Starting from the first telescopic sunspot observations by David and Johannes Fabricius, 
Galileo Galilei, Thomas Harriot and Christoph Scheiner, the 400-year sunspot record is 
one of the longest directly recorded scientific data series, and forms the basis for numerous 
studies in a wide range of research such as, e.g., solar and stellar physics, solar-terrestrial 
relations, geophysics, and climatology. During the 400- year interval, sunspots depict a grea t 



deal of variability from the extre mely quiet period o f the Maunder minimum (lEddy 



to the very active modern time (jSolanki et al 



1976|) 



20041 ). The sunspot numbers also form a 



benchmark data series, upon which virtually all modern models of long-term solar dynamo 
evolution, either theoretical or (semi)empirical, are based. Accordingly it is important to 
review the reliability of this series, especially since it contains essential uncertainties in the 
earlier part. 

The first sunspot number series was introduced by Rudolf Wolf who observed sunspots 
from 1848 until 1893, and const ructed the monthly sunspot numbers since 1749 using 



archival records and proxy data (jWoli 



18611 ). Sunspot activity is dominated by the 11-year 
cyclicity, and the cycles are numbered in Wolf's series to start with cycle #1 in 1755. 
When constructing his sunspot series Wolf interpolated over periods of sparse or missing 



sunspot observations u sing geomagnetic 



evolution of sunspots (jHoyt et al 



1994 



proxy data, thus 


osing 


Hovt & Schatten 


1998) 



osing the actual detailed temporal 



19981 ). Sunspot observations were 



particularly sparse in the 1790's, during solar cycle #4 which became the longest solar 
cycle in Wolf's reconstruction with an abnormally long declining phase (see Fig. [T]A.). The 
quality of Wolf's sunspot series during that period has been questioned since long. Based 
on independent auroral observations, it was proposed by Elias Loomis already in 1870 that 
one sm all solar cycle may have been completely lost in Wolf's sunspot reconstruction in the 



1790's ( Loomis 



18701 ). being hidden inside the interpolated, exceptionally long declining 
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phase of solar cycle #4. This extraordinary idea was not accepted at that time. A century 
later, possible errors in Wolf's compilation for the late 18 th century have been emphasize d 



again based on detailed studies of Wolf's sunspot series (jGnevyshev &: Ohl 



1983|). 



1948 



Sonett 



Recently, a more extensive and consiste nt sunspot number series (Fig. [TJB), the group 



sunspot numbers (GSN), was introduced by 



Hoyt Schattenl (119981 ) . which increases 



temporal resolution and allows to evaluate the statistical uncertainty of sunspot numbers. 
We note that the GSN series is based on a more extensive database than Wolf's series and 
explicitly includes all the data collected by Wolf. However, it still depicts large data gaps 
in 1792-1794 (t his interval wa s interpolated in Wolf's series). Based on a detailed study of 



the GSN series. 



Usoskin et al. 



(12001b 



2002) revived Loomis' idea by showing that the lost 



cycle (a new small cycle started in 1793, which was lost in the conventional Wolf sunspot 
series) agrees with both the GSN data (Fig. [T]B) and indirect solar proxies (aurorae) and 
does not contradi ct with the cosmogen ic isotope data. The existence of the lost cycle has 



been disputed by iKrivova et al. 



numbers. However, gued by 



(200 21) based on data of cosm ogenic isotope and sunspot 



Usoskin et al. 



(12002 . 



20031 ). the lost cycle hypothesis 



does not contradict with sunspots or cosmogenic isotopes and is supported by aurorae 
observations. Using time series analysis of sparse sunspot counts or sunspot proxies, it is 
hardly possible to finally verify the existence of the lost solar cycle. Therefore, the presence 
of the lost cycle has so far remained as an unresolved issue. 

Here we analyze newly restored original solar drawings of the late 18th century to 
ultimately resolve the old mystery and to finally confirm the existence of the lost cycle. 
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2. Data and analysis 



2.1. Positional sunspot data 



Most of Wolf's sunspot numbers in 1749-1796 were constructed from observations by 
the German amateur astronomer Johann Staudacher who not only counted sunspots but 
also drew solar images in the second half of the 18th century (see an example in Fig. [2]). 
However, only sunspot counts have so far been used in the sunspot series, but the spatial 
distribution of spots in these drawings has not been analyzed earlier. The first analy sis o 



this data, which covers the lost cycle period in 1790's, has been made only recently (Arlt 
20081 ) using Staudacher's original drawings. Additionally, a few original solar disc drawings 
made by the Irish astronomer James Archibald Hamilton and his assistant s ince 1795 have 



2009bh. After t 



lie 



2008 



2009a 



b) 



been recently found in the archive of the Armagh Observatory ( lArlt 
digitization and processing of these two sets of original drawings (lArlt 
the location of individual sunspots on the solar disc in the late 18th century has been 
determined. This makes it possible to construct the sunspot butterfly diagram for solar 
cycles #3 and 4 (Fig. [EC), which allows us to study the existence of the lost cycle more 
reliably than based on sunspot counts only. 



2.2. Constructing the butterfly diagram from data with uncertainties 



Despite the good quality of original d rawings, ther e is an uncertainty in determining 



the actual latitude for some sunspots (see 



Arlt 



(12009a 



tj) for details). This is related to 



the limited information on the solar equator in these drawings. The drawings which are 
mirrored images of the actual solar disc as observed from Earth, cannot be analyzed by 
an automatic prodecure adding the heliographic grid. Therefore, special efforts have been 
made to determine the solar equator and to place the grid of true solar coordinates for 
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each drawing (see Fig. [2]). Depending on the information available for each drawing, the 
unce rtainty in d efining the solar equator, AQ, ranges from almost 0° up to a maximum of 



15° f lArlt 



2009al ). The latitude error of a sunspot, identified to appear at latitude A, can be 
defined as 

A A = AQ • sin(a) (1) 

where AQ is the angular uncertainty of the solar equator in the respective drawing, and a 
is the angular distance between the spot and the solar disc center. Accordingly, the final 
uncertainty A\ can range from 0° (precise definition of the equator or central location of 
the spot) up to 15°. We take the uncertain spot location into account when constructing 
the semiannually averaged butterfly diagram as follows. Let us illustrate the diagram 
construction for the second half-year (Jul-Dec) of 1793 (Fig. [3]). During this period there 
were only two daily drawings by Staudacher with the total of 8 sunspots: two spots on 
August 6th, which were located close to the limb near the equator, and six spots on 
November 3rd, located near the disc center at higher latitudes. The uncertainty in definition 
of the equator was large (AQ = 15°) for both drawings. Because of the near-limb location 
(large a) of the first two spots, the error A a of latitude definition (Eq. [1]) is quite large. 
The high-latitude spots of the second drawing are more precisely determined because of 
the central location of the spots. The latitudinal occurrence of these eight spots and their 
uncertainties are shown in Fig. [3] as stars with error bars. The true position of a spot is 
within the latitudinal band A ± A\, where A\ is regarded as an observational error and A 
as the formal center of the latitudinal band. Accordingly, when constructing the butterfly 
diagram, we spread the occurrence of each spot within this latitudinal band with equal 
probability (the use of other distribution does not affect the result). Finally, the density 
of the latitudinal distribution of spots during the analyzed period is computed as shown 
by the histogram in Fig. [3j This density is the average number of sunspots occurring per 
half-year per 2° latitudinal bin. Each vertical column in the final butterfly diagram shown 
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in Fig. [Ip is in fact such a histogram for the corresponding half-year. 



2.3. Statistical test 



Typically, the sunspots of a new cycle appear at rather high latitudes of about 
20-30°. This takes place around the solar cycle minimum. Later, as the new cycle evolves 
the sunspot emergence zone slowly moves towards the solar equator. Thi s recurren t 



"butterfly" -like pattern of sunspot occurrence is known as t 
and is related to the action of the solar dynamo (see, e.g. 



re Sporer law 



Charbonneau 



Maunder 



190J) 



20051 ). It is 



important that the systematic appeara nce of sunspots a t high latitudes unambiguously 



19751 ) and thus may clearly distinguish 



indicates the beginning of a new cycle ( IWaldmeier 
between the cycles. 

One can see from the reconstructed butterfly diagram (Fig. [TP) that the sunspots in 
1793-1796 appeared dominantly at high latitudes, clearly higher than the previous sunspots 
that belong to the late declining phase of the ending solar cycle #4. Thus, a new "butterfly" 
wing starts in late 1793, indicating the beginning of the lost cycle. 

Since sunspot observations are quite sparse during that period, we have performed a 
thorough statistical test as follows. The location information of sunspot occurrence on the 
original drawings during 1793-1796 (summarized in Tabled]) allows us to test the existence 
of the lost cycle. The observed sunspot latitudes were binned into three categories: low 
(< 8°), mid- (8°— 16°) and high latitudes (> 16°), as summarized in column 2 of Table[2j We 
use all available data on latitude distribution of sunspots since 1874 covering solar cycles 12 
through 23 (The combined Royal Greenwich Observatory (1874-1981) and USAF/NOAA 



(1981-2007) Sunspot Data set: http://solarscience.msfc.nasa.gov/greenwch.shtml) as the 
reference data set. We tested first if the observed latitude distribution of sunspots (three 
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daily observations with low-latitude spots, one with mid-latitude and three with only 
high-latitude spots, see Table [2]) is consistent with a late declining phase (D-scenario, 
i.e. the period 1793-1796 corresponds to the extended declining phase of cycle #4) or 
with the early ascending phase (A-scenario, i.e., the period 1793-1796 corresponds to the 
ascending phase of the lost cycle). We have selected two subsets from the reference data set: 
D-subset corresponding to the declining phase which covers three last years of solar cycles 
12 through 23 and includes in total 11235 days when 33803 sunspot regions were observed; 
and A-subset corresponding to the early ascending phase which covers 3 first years of solar 
cycles 13 through 23 and includes 10433 days when 47096 regions were observed. 

First we analyzed the probability to observe sunspot activity of each category on a 
randomly chosen day. For example, we found in the D-subset 4290 days when sunspots were 
observed at low latitudes below 8°. This gives the probability p = 0.38 (see first line, column 
3 in Table [2]) to observe such a pattern on a random day in the late decline phase of a cycle. 
Similar probabilities for the other categories in Table [2] have been computed in the same 
way. Next we tested whether the observed low-latitude spot occurrence (three out of seven 
daily observations) corresponds to declining/ascending phase scenario. The corresponding 
probability to observe n events (low-latitude spots) during m trials (observational days) is 
given as 

p m =p n -a-p) m - n -c m , (2) 

where p is the probability to observe the event at a single trial, and is the number of 
possible combinations. We assume here that the results of individual trials are independent 
on each other, which is justified by the long separation between observational days. Thus, 
the probability to observed three low-latitude spots during seven random days is = 0.27 
and 0.07 for D- and A-hypotheses, respectively. The corresponding probabilities are given 
in the first row, columns 5-6 of Table [2J The occurrence of three days with low latitude 
activity is quite probable for both declining and ascending phases. Thus, this criterion 



- 9- 



cannot distinguish between the two cases. The observed mid-latitude spot occurrence (one 
out of seven daily observations) is also consistent with both D- and A-scenarios. The 
corresponding confidence levels (0.06 and 0.22, respectively, see the second row, columns 5-6 
of Table [2]) do not allow to select between the two hypotheses. Next we tested the observed 
high-latitude spot occurrence (three out of seven observations) in the D/ A-scenarios (the 
corresponding probabilities are given in the third row of Table [2]). The occurrence of 
three days with high-latitude activity is highly improbable during a late declining phase 
(D-scenario). Thus, the hypothesis of the extended cycle #4 is rejected at the level of 
5 • 10~ 4 . The A-scenario is well consistent (confidence 0.26) with the data. Thus, the 
observed high-latitude sunspot occurrence clearly confirms the existence of the lost cycle. 

We also noticed that sunspots tend to appear in Northern hemisphere (13 out of 
16 observed sunspots appeared in the Northern hemisphere). Despite the rather small 
number of observations, the statistical significance of asymmetry is quite good (confidence 
level 99%), i.e. it can be obtained by chance with the probability of only 0.01, in a 
purely symmetric distribution. Nevertheless, more data are needed to clearly evaluate the 
asymmetry. 

Thus, a statistical test of the sunspot occurrence during 1793-1796 confirms that: 

• The sunspot occurrence in 1793-1796 contradicts with a typical latitudinal pattern in 
the late declining phase of a normal solar cycle (at the significance level of 5 • 10~ 4 ). 

• The sunspot occurrence in 1793-1796 is consistent with a typical ascending phase 
of the solar cycle, c onfirming the start of the lost solar cycle. We note that it has 



been shown earlier (jUsoskin et al. 



20031 ). using the group sunspot number, that the 
sunspot number distribution during 1792-1793 was statistically similar to that in the 
minimum years of a normal solar cycle, but significantly different from that in the 
declining phase. 
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• The observed asymmetric occurrence of sunspots during the lost cycle is statistically 
significant (at the significance level of 0.01). 



Therefore, the sunspot butterfly diagram (Fig. [IP) unambiguously proves the existence 



of the los t cycle in the 



numb ers ( lUsoskin et al. 
2002k 



ate 18 



; h cent ury, verifying the earli er evidence based on sunspot 



2001b . 



20031 ) and aurorae borealis (ILoomis 



1870 



Usoskin et al. 



3. Discussion and Conclusions 



An additional cycle in the 1790's changes cycle numbering before the Dalton 
minimum, thus verifying the validity of the Gnevyshev-Ohl rul e of sunspot cycle p airing 



(iGnevyshev fc Ohl 



20011) in 



19481 ) and the related 22-year periodicity (jMursula et al. 
sunspot activity throughout the whole 400-year interval. Another important consequence 
of the lost cycle is that, i nstead of one abnor mally long cycle #4 (min-to-min length ^15.5 



years according to GSN (lUsoskin et al. 



20021 ) ) there are two shorter cycles of about 9 and 7 



years (see Fig. HP). Note also that some physical dynamo models even p redict the existenc e 



19991). 



of cycles of small amplitude and short duration near a grand minimum (jKiiker et al. 
The cycle #4 (1784-1799 in GSN) with its abnormally long duration dominates empirical 
studies of relations, e.g., between cycle length and amplitude. Replacing an abnormally long 
cycle #4 by one fairly typical and one small short cycle changes empirical relations based on 
cycle length statistics. Thi s will affect , e.g., predictions of f uture solar activity by statistical 



or dynamo-based models (IDicke 



1978 



Dik pati et al 



2006 



Bra isa et al 



impo rtant solar-terrestrial relations (IFriis-Christensen fc Lassen 
1992k 



1991 



2009 ), and som e 



Kellv & Wiriev 



The lost cycle starting in 1793 depicts notable hemispheric asymmetry with most 
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sunspots of the new cycle occurring in the northern solar hemisphere (Fig. (DC). This 
asymmetry is statistically significant at the confidence level of 99%. A similar, highly 
asymmetric sunspot distribution existe d during the Maunder minimum of sunspot activi ty 



in the second half of the 17th century (ISoon fc Yaskel 



2003 



Ribes fc Nesme-Ribes 



19931 ). 



However, the suns pots during the Maunder min imum occurred preferably in the southern 



solar hemisphere (ISokoloff fc Nesme-Ribes 



1994 ). i.e., opposite to the asymmetry of the 



lost cycle. This shows that the asymme try is not constant, con t rary to some earlier models 



involv ing the fossil solar magnetic field (jBravo fc Stewart 



1995 



Boruta 



1996 



Usoskin et al. 



2001al ). Interestingly, this change in hemispheric asymmetry between the Maunder and 
Dalton minimum is in agreement with an earlier, independent observation, based on 
long-term geomagnet ic activity, that the nort h-south asymmetry oscillates at the period of 



about 200-250 years (IMursula fc Zieger 



200l|). 



Concluding, the newly recovered spatial distribution of sunspots of the late 18th 
century conclusively confirms the existence of a new solar cycle in 1793-1800, which has 
been lost under the preceding, abnormally long cycle compiled by Rudolf Wolf when 
interpolating over the sparse sunspot observations of the late 1790's. This letter brings 
the attention of the scientific community to the need of revising the sunspot series in the 
18th century and the solar cycle statistics. This emphasizes the need to search for new, 
yet unrecovere d, solar data to restore details of solar activity evolution in the past (e.g., 



Vaquero 



20071 ). The new cycle revises the long-held sunspot number series, restoring its 



cyclic evolution in the 18th century and modifying the statistics of all solar cycle related 
parameters. The northern dominance of sunspot activity during the lost cycle suggests 
that hemispheric asymmetry is typical during grand minima of solar activity, and gives 
independent support for a systematic, century-scale oscillating pattern of solar hemispheric 
asymmetry. These results have immediate practical and theoretical consequences, e.g., to 
predicting future solar activity and understanding the action of the solar dynamo. 
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Table 1: Sunspot occurrence during the lost cycle with dates, number (N) and latitude 



range (A) of the observed spots on each day. 



Date and Observer 


N 


A 


1793.08.06 Staudacher 


2 


< 3° 


1793.11.03 Staudacher 


6 


> 18° N 


1795.02.19 Staudacher 


3 


> 20° N 


1795.10.08-15 Hamilton 


2 


15° S and 6° N 


1795.11.02-03 Hamilton 


1 


5.5° S 


1796.01. 31 a Staudacher 


2 


> 16° N 



"Shown in Fig. H 



Table 2: Number of days {n) with observed sunspots in 1793-1796, sunspot latitude ranges, 
probability (p) of sunspots to be found on a randomly selected day, and the cumulative 
probability (P), calculated for D- and A-scenaria using all data since 1874. 



n 


Latitude range 


probability p 
D-scenario A-scenario 


cumulative 
D-scenario 


probability P 
A-scenario 


3 


Low latitude (< 8°) 


0.38 


0.1 


0.27 


0.07 


1 


Mid-latitude (8-16°) 


0.48 


0.32 


0.06 


0.22 


3 


High latitude (> 16°) 


0.026 


0.52 


0.0005 


0.26 
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Fig. 1. — Sunspot activity in late 18th and early 19th century. The start of the lost cycle 
(late 1793) is denoted by a vertical dashed line. A: Monthly Wolf sunspot numbers with 
conventional solar cycle numbers shown on the top. B: Monthly group sunspot numbers. 
C: The newly reconstructed sunspot butterfly diagram, which takes into account the uncer- 
tainties of the estimate sunspot latitudes. The color scale on the right gives the density (in 
year~ 1 deg~ 1 ) of sunspots in latitude-time bins (one bin covers 2° in latitude and 6 months 
in time). Grey bars indicate that no latitudinal information is available for the correspond- 
ing half-year. D: Lengths of solar cycles. The conventional lengths using the group sunspot 



Fig. 2. — An example of drawings of sunspots on the solar disc made by Staudacher in 
Niirnberg, Germany at 1 p.m. local time on January 31, 1796. The image is mirrored to 
correspond to the real view to the solar disc. The heliographical grid has been included 



during the image processing (lArlt 



2009ah . 
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Fig. 3. — An example of the sunspot latitudinal distribution for Jul-Dec 1793, with 2° 
latitudinal bins. Stars with error bars denote latitudes of individual spots as defined from 
Staudacher's drawings. The histogram depicts the density of the latitudinal distribution of 
sunspots per 2° bins per half-year computed by including the uncertainties. 



